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ABSTRACT 

Since 2005, the blazar 3C 454.3 has shown remarkable flaring activity at all 
frequencies, and during the last four years it has exhibited more than one 7-ray 
flare per year, becoming the most active 7-ray blazar in the sky. We present 
for the first time the multi- wavelength AGILE, Swift, INTEGRAL, and GASP- 
WEBT data collected in order to explain the extraordinary 7-ray flare of 3C 454.3 
which occurred in November 2010. On 2010 November 20 (MJD 55520), 3C 454.3 
reached a peak flux (E>100 MeV) of = (6.8±1.0) x 10"^ photons cm'^ g-^ on 
a time scale of about 12 hours, more than a factor of 6 higher than the flux of the 
brightest steady 7-ray source, the Vela pulsar, and more than a factor of 3 brighter 
than its previous super-flare on 2009 December 2-3. The multi-wavelength data 
make a thorough study of the present event possible: the comparison with the 
previous outbursts indicates a close similarity to the one that occurred in 2009. 
By comparing the broadband emission before, during, and after the 7-ray flare, 
we find that the radio, optical and X-ray emission varies within a factor 2-3, 
whereas the 7-ray flux by a factor of 10. This remarkable behavior is modeled 
by an external Compton component driven by a substantial local enhancement 
of soft seed photons. 

Subject headings: galaxies: active - galaxies: quasars: general - galaxies: quasars: 
individual: 3C 454.3 - galaxies: jets - radiation mechanism: non thermal 
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Introduction 



The flat spectrum radio quasar 3C 454.3 (PKS 2251+158; z = 0.859) is the brig htest 
7-ray ( 0.1-10 GeV) blazar detected after the launch of the AGILE (iTavani et al.ll2009l ) and 
Fermi flAtwood et al.lboOQl) satelhtes. During 2007-2010, AGILE detected and investigate d 
several gamma-ray flares ( iVercellone et al.ll2010al : iPacciani et al.ll2010l : IStriani et al.lbOlOdl ). 
These observations allowed us to establish a possible correlation between the 7-ray (0.1 -10 
GeV) and the optical (R band) flux variations with no time delay, or with a lag of the former 
with respect to the latter of about half a day. Moreover, the detailed physical modeling of 
the spectral energy distributions (SEDs) when 3C 454.3 was at different flux levels provided 
an interpretation of the emission mechanism responsible for the radiation emitted in the 7- 
ray energy band, assumed to be inverse Compton scattering of photons from the broad line 
region (BLR) clouds off the relativistic electrons in the jet, with bulk Lorentz factor F ~ 20. 
Similar re sults where obtained by other groups, by analyzing Fermi and multi-wavelength 



data (e.g..lGhisellini et al.ll2007l : iBonning et al.ll2009l : lAbdo et al.ll2009l : iFoschini et al.ll2010 



Bonnoli et al.ll201ll ). During the period 2009 December 2-3, 3C 454.3 exhibited an intense 
7-ray flare reaching a peak value of _pp>2009 = (2.0 ±0.4) x 10^^ photons cm^^ s^^ in the range 
0.1-10 GeV. This extreme behavior in the 7-ray might require a more sophisticated modeling 
with respect to the widely a ccepted one-zone, synchro tron self-Compton (SSC) and exter- 
nal Compton (EC) models (jCelotti fc Ghisellinil |2008| . for a review of the blazar emission 
mechanisms and energetics). In particular, the almost simultaneous multi- wavelength data 
collected during the 7-ray flare required an additional popul ation of accelerated e lectrons 
co-existing with the soft seed photons of the SSC/EC model (IPacciani et al.ll2010l ). Alter- 
natively to the S SC/EC models, the blazar SEP can be explained in the framework of t he 
hadronic models (IMannheim &: Schlickeiserlll994j : iMiicke &: Protherodl200ll : iBeckerl 120081 ) . 



Recently, iFoschini et al.l (120 111 ) and lAbdo et al.l (l201ll ). analyzing Fermi data acquired 
during the November 2010 flare, show that the extremely fast variability, on a time-scale of 
about 3-6 hr, favors a 7-ray emission originating from a compact region, below the pc scale. 

In this Letter we report on the multi- wavelength AGILE, Swift, INTEGRAL, and GASP- 
WEBT campaign covering the extraordinary November 2010 7-ray flare of 3C 454.3. The 
quoted uncertainties are given at the la level, unless otherwise stated, and a A-CDM cos- 
mology {h = 0.71, Qrn = 0.27, and = 0.73) was adopted. 
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2. Data Analysis 



AGILE detected increased and prolonged 7-ray emission from 3C 454.3 starting from 
2010 October 28 (MJD 55497), with a maxim um emission on 2010 November 20 (MJD 55520, 
Vercellone et al.lboiObi Istriani et aPboiOal lblD). L evel-1 AGILE-GRID da ta were analyzed 
using the AGILE Standard Analysis Pipeline (see IVercellone et al.ll2010ai for a description 
of the AGILE data reduction). We used 7-ray events filtered with the FM3 . 119 AGILE Filter 
pipeline. Counts, exposure, and Galactic background 7-ray maps were created with a bin-size 
of 0?5 X 0?5 , E > 100 MeV, and including all events collected up to 60° off-axis. To reduce 
the 7-ray Earth albedo contamination, we rejected all 7-ray events whose reconstructed 
directions form angles with the satellite-Earth vector < 85°. We used the latest version 
(BUILD-20) of the Calibration files (lOOlO), which will be pubhcly availa ble at the ASI 
Scien ce Data Centre (ASDC) sit^, and of the 7-ray diffuse emission model ( Giuliani et al. 



2OO4J ). We ran the AGILE Multi-Source Maximum Likelihood Analysis (ALIKE) task with 
an analysis radius of 10°. 

The Swift X-ray Telescope (XRT) data were processed with standard procedures (xrt- 
PIPELINE vO.12.6 within HEASOFT V.6.10), filtering and screening criteria. The XRT data 
were in windowed timing mode (WT). Spectra were extracted on an orbit-by-orbit basis due 
to a rotation of the roll angle, within circular regions with a radius of 47.15 arcsec. We 
used the latest spectral redistribution matrices (20101206). The Swift Ultraviolet and Op- 
tical Telescope (UVOT) data analysis was performed using the uvotimsum and uvotsource 
tasks. Source counts were extracted from a circular region with a 5 arcsec radius. The 
background was extracted from a nearby source-f ree circular r egiori s. The reported fluxes 
are on the UVOT photometric system described in lPoole et al.l (120081 ) . and are not corrected 
for Galactic extinction. 

The INTEGRAL data were obtained as a Target of Opportu nity (ToO) observ ation trig- 
gered immediately after the 7-ray flare (for preliminary results see lPian et al.ll2010r). and were 



proce ssed using the ISDC Off-line Scientific Analysis software (OSA V . 9 . 0, ICourvoisier et al. 



20031 ). INTEGRAL observed the source between 2010 November 21.70 UT and November 
27.24 UT, for a total on source time of 400 ks. Light curves and spectra were extracted for 
each individual science window (SCW) and later combined. We investigated the 4-26 keV 
(JEM-X both units) and 20 - 200 keV (IBIS/ISGRI) energy ranges. We extracted the hght- 
curve from IBIS/ISGRI in the 20-40 keV, and 40-100 keV energy bands to verify possible 
variability on time scales from hours to one day: the average count rate in the softer (harder) 
band is 1.6±0.1 counts s~^ (1.4 ±0.1 counts s~^), or about 11 (15) mCrab without significant 



^http : //agile . asdc . asi . it| 
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Fig. 1. — Multi- wavelength light curves from radio to 7-ray. 7-ray (E>100MeV) data are 
in units of 10~^ photons cm~^ s~^; hard X-ray (20-100 keV) data are in units of counts s~^; 
soft X-ray (2-10 keV) data are in units of 10~^^ erg cm~^ s~^; UV, optical and NIR data 
are in units of mJy; radio data are in units of Jy. See §[3] for details. The starting date, 
MJD 55500, corresponds to 2010-10-31 00:00 UT. 
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variability. We then analyzed the averaged spectrum using 8 pre-defined energy channels in 
JEM-X and 16 customized energy channels for IBIS/ISGRI with XsPEC (v. 12.6.0). We as- 
sumed a single power law model and free cross-calibration constants for JEM-X (Ci and C2) 
with respect to ISGRI, and a 3% systematic error to account for the instrumental calibration 
uncertainties. We obtain (uncertainties at 90% c.l.): T = 1.59 ±0.08 {Xred/^-'^-^- — 1-35/23), 
Ci = 1.1 ± 0.2 , C2 = 1.3 d= 0.2 , F4„200kcv = (4.7t°-2^ 



-oi) X 10- 



1 n —2 
erg cm s 



The GLAST- AGILE Support Program (GASP; IVillata et al.ll2008l 120091 ) is a project 
born from the Whole Earth Blazar Telescop^ (WEBT) in 2007, which has analys ed the 



multi - frequ e ncy behavior of 3C 4 54.3 sin ce the unprecedented optical outburst of 2 005 (IVillata et al 



2006. 


2007; 


Raiteri et al. 


2007. 


2008aUbl: 


Villata et al. 


2009; 


Raiteri et al. 


2011) 



band GASP observations of 3C 454.3 in the period considered in this paper were performed 
by the following observatories: Abastumani, Calar Alto, Crimean, Galaxy View, Goddard 
(GRT), Lowell (Perkins), Maria Mitchell, New Mexico Skies, ROVOR, Roque de los Mucha- 
chos (KVA an d Liverpool) , Saba dell, and St. Petersburg. Additional da ta were taken by 
the A AVSO§ JKraici et al.lboiol l and the Yale Fermi/SMARTS projecfl dchatteriee et al. 



2OIII ). Some of the above observatories also provided data in the B, V, and / bands. Near- 



IR (NIR) data in J, H, and K bands are from Campo Imperatore. Radio flux densitie s 
were measured at: Submillimeter Array (SMA, 345 and 230 GHz, see iGurwell et al.l 120071 ). 
Medic ina (8 GHz) and UM RAO (14.5, 8.0, and 4.8 GHz). Data reduction and analysis 
follow iRaiteri et all f|2008al ) . 



3. Results 

Figure [T] shows the multi-wavelength light curves from radio to 7-ray; from top to bot- 
tom: AGILE/GRID (E>100 MeV) at a variable time-bin of^ 2, 1, and 0.5-day; INTEGRAL 
(20-100 keV); Swift/XRT (2-10 keV); Swift/UVOT UV wl (red circles), m2 (green squares), 
and w2 (blue triangles); Swift/UVOT (triangles), and GASP- WEBT (squares) optical U 
(blue points), B (green points), and V (red points); GASP- WEBT optical R (blue triangles, 
with additional data from AAVSO as green squares) and I (red circles); GASP- WEBT near 
infra-red J (red points), H (green points), and K (blues points); and GASP- WEBT 4.8 GHz 
(red circles), 8 GHz (green circles), 14.5 GHz (blue circles), 230 GHz (black circles), and 
345 GHz (yellow circles). 



^http : //www . oato . inaf . it/blazairs/webt/ 
" http : //www . aavso . org/| 

' http : //www . astro . yale . edu/ smarts/ glast/ 
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Fig. 2. — Zoom of the 7-ray (top panel), i?-band (middle panel), and 230 GHz (bottom 
panel) light curves. Units and energy bands are the same as in Figure [H The starting date, 
MJD 55515, corresponds to 2010-11-15 00:00 UT. 
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When data are available, a clear peak is present at approximately MJD 55520 in all 
light curves. In particular, the 7-ray light curve shows different behavior during the pre- and 
the post-ffare periods, the former being much steadier than the latter, on almost the same 
time-scale. The UV and optical light curves show a very similar trend, in particular the 
remarkably fast flare centered at approximately MJD 55510, with rising and falling of about 
a factor of 2-2.5 in about 48 hours. During the same period, the X-ray flux varied by about 
20%, while no significant variability is detected at other wavelengths. We also note that 
there is an average offset between the GASP-WEBT and Swift/UVOT data in the V and B 
band. If (V, B) and (f , h) are the GASP-WEBT and Swift/UVOT data , respectively, then a 
good agreement is found when B — b = 0.1 mag, and V — v = —0.05 mag. iRaiteri et al.l (1201 ll ) 
discuss the nature of this offset and provide a tool for the calibration of the Swift/UVOT 
data. Figure [2] shows a zoom, centered on the 7-ray flare date, of the 7-ray, i?-band and 
230 GHz light curves. We note that, to the sudden enhancement (about a factor of 4) of the 
gamma-ray flux on MJD 55517, does not correspond a similar flare in the other wavebands. 

The Swift/yiKT spectra were rebinned to have at least 20 counts per en ergy bin, and 



were fit with an absorbed power law model. Following I Vercellone et al.l (l2010al) . the Galactic 
absorption was fixed to the value of A^^^' = 1.34 x 10^^ cm~^ ( Villata et al. I2OO6 ). In Figure[2] 



(a), we show the X-ray photon index as a function of the 2-10 k eV X-ray flux. N o partic ular 
trend is present, contrary to what was previously reported in I Vercellone et al.l (l2010al ) for 
the period 2007-2009. 

The 7-ray spectra were computed in four different time periods, pre-flare (MJD 55497.60 
- 55516.40), flare (MJD 55518.25 - 55520.25), post-flare (MJD 55521.25 - 55543.75) and 
whole period. Table [1] shows the spectral parameters that we obtained by fitting the AGILE 
data with a simple power-law: 



F{E) = kx 



IMeV J 



(1) 



Table H] reports the numerical values during the different time periods. We note that the flux 
during the flare is about a factor of 2 5 larger than the average value reported in the First 
AGILE Catalogue jPittori et al.lbood . July 2007 - June 2008). 



4. Discussion 



Since its launch in 2007, AGILE detected significant 7-ray emission from 3C 454.3 with 
repeated and prolonged flaring activity. Figure |3](b), shows the ratio -^highZ-^iow ^ function 
of time, where F^-^ is the flux maximum value during each flare detected by AGILE, and 
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Fig. 3. — Panel (a): X-ray photon index as a function of the X-ray flux in the energy band 
2-10 keV. Panel (b): peak over low 7-ray flux ratio for the major AGILE detected flares as 
a function of time. Panel (c): blue circles represent the ratio between the 7-ray and i?-band 
peak fluxes as a function of the i?-band peak flux, while red square represents the same 
quantity for the fast optical-UV flare on MJD 55510. 



- 11 - 



-F[^^ is the lowest flux point just before the corresponding rising in the light curve. While in 
2007-2008 this ratio is of the order of 3-5, in 2009 and 2010 this ratio is at least a factor 
of 2-3 larger. It is worth noting that, during this latest flare, the 230 GHz, optical-UV, 
X-ray, and 7-ray flux variations are almost simultaneous. The 230 GHz flux shows a more 
prolonged active state after the 7-ray super-flare (~ 10 days), while the UV and optical 
fluxes reach levels comparable or lower than the pre-flare ones. The extremely large 7-ray 
dynamic range in the 2010 flare, the high 230 GHz flux level, and the relatively bright state 
in the optical band, when compared with other 7-ray flares, could not be explained only i n 



terms of the alignment of different regions of the jet as suggested in IVercellone et al.l (l2010al ). 
Figure [3] (c), provides further evidence of the different behavior of the 2009 and 2010 flares 
with respect to the previous ones. The plot of the ratio between the 7-ray and -R-band peak 
fluxes as a function of the i?-band peak flux (blue circles) shows a difference among flares in 
the different epochs. 

The red square represents the point relative to the fast optical-UV flare on MJD 55510. 
This last point is particularly interesting because of the fast rising and decay time scales 
(~ld). To explain the lack of a simultaneous 7-ray flare, a magnetic field enhancement by 
a factor of 2, a local 7-7 absorption, or a locally seed-photon starved zone are required. 

The INTEGRAL 20-100 keV light curve shows another peculiarity of the November 2010 
flare with respect to the intense optical and X-ray flare of May 2005. During the November 
2010 flare, the source shows only moderate flux variability (~ 1.6a level) and a relatively 
hard spectrum (F = 1.59 ± 0.08), while in 2005, during a comparably long observation at a 
comparable flu x level, the source exhibited remarkable flux variability and a softer spectrum 



(F = 1.8 ± 0.1, iPian et al.ll2006l ). Despite the moderate variability in the 20-100 keV band. 
Figure [1] seems to show a similar trend between the light curves in the soft, hard X-ray and 
7-ray bands. 

Figure H] shows SEDs chosen at three epochs with the largest number of data points. 
The NIR, optical, UV, and X-ray data are corrected for Galactic extinction. Green open 
triangles represent the pre-flare SED, accumulated on MJD 55512. In order to obtain enough 
statistics, the AGILE data were acquired in the period MJD 55497.60 - 55516.40. Red filled 
circles represent the flare SED, accumulated on MJD 55519 (AGILE data: MJD 55518.25 
- 55520.25). Blue open squares represent the post-flare SED, accumulated on MJD 55526 
(AGILE data: 55521.25 - 55543.75). The INTEGRAL IBIS/ISGRI integrated spectrum 
(MJD 55521.70 - 55527.30) is also reported. 

We fit the pre-flare SED by means of a one- zone leptonic model, considering the contri- 
butions from SSC and from external seed photons originating both fr om the accretion dis k 



(AD) and from the BLR (detailed description of this model is given in lVittorini et al.ll2009l ). 



Table 1. Top: 7-ray spectral parameters. Bottom: parameters for the pre-flare and flare, 

SED models, respectively (see §|l]for details). 



Parameter Pre-flare Flare Post-flare Whole period Units 



7-ray spectral parameters 


Fe>100 MeV 


6.2 ±0.6 


53.5 ±3.7 


13.8 ±0.6 16.^ 


^±0.7 


10 ^phcm ^s ^ 


k 


0.4 


11.1 


10.0 


3.8 


10-3 ph cm-^s-^MeV"^ 


a 


1.93 ±0.20 


2.13 ±0.13 


2.37 ±0.08 2.15 ±0.08 








SEDs 


model parameters 






ai 


2.35 


2.35 








ah 


4.2 


4.8 








7min 


50 


80 








7b 


650 


700 








K 


300 


700 








Rblob 


7.0 


3.6 






10^^ cm 


B 


0.65 


1.1 






G 
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Indeed, emission from bot h the BLR and the AD were detected during faint states of the 
source (iRaiteri et al.l 120071 ). The sohd hues represent the total contributions before (green), 
during (red), and after the flare (blue). 



Hadronic models are chal 
wavelength, as p ointed out by 



AGILE (but see lBarkov et al.l 120101 ). 



enged by the de t ection of correlated variability at different 



Finke fc Dermed (120101 ). and as observed by both Fermi and 



The emission along the jet is assumed to be produced in a spherical blob with comoving 
radius -Rbiob by accelerated electrons characterized by a comoving broken power law energy 
density distribution. 



^e(7) 



(2) 



(7/7b)"' + il/lh)"^ ' 

where 7 is the electron Lorentz factor varying between 80 < 7 < 8 x 10^, ai and are the 
pre- and post-break electron distribution spectral indices, respectively, and 7b is the break 
energy Lorentz factor. We assume that the blob contains an homogeneous and random 
magnetic field B and that it moves with a bulk Lorentz Factor F at an angle Go with respect 
to the line of sight. The relativistic Doppler factor is 5 = [F (1 — /3 cosGq)]"^, where (3 is 
the blob bulk speed in units of the speed of light. Our fit parameters are listed in Table [TJ 

Our modeling of the 3C 454.3 high-energy emission is based on an inverse Compton 
(IC) model with two main sources of external seed photons: (1) the AD characterized by a 
blackbody spectrum peaking in the UV with a bolometric luminosity for an IC-scattering 
blob at a distance from the central part of the dislj^; (2) the BLR with a spectrum 



peaking in the V, placed at a distance from the central black ho 
and assumed to reprocess 10% of the irradiating co ntinuum (see 



e of rRT;R = 3 X 10 



18 



cm. 



Vercellone et al. 2009 for 



2010 for the December 



the first application of this model to 3C 454.3, and iPacciani et al. 
2009 fiare modeling). 

The fiaring behavior in the optical and 7-ray energy bands is puzzling. We observe a 
first rapid fiare at MJD 55510 with the optical fiux rising and falling by a factor 2 in about 
24 hours without a 7-ray counterpart. After a time T=7d, at MJD 55517, the optical fiux 
increases by a factor of 2, followed after about 12 hours by a 7-ray increase by a factor of 4. 
On the contrary, on MJD 55520, the fiux variations in the optical and 7-ray bands are of the 
same magnitude, as expected from an EC mechanism. This complex behavior challenges the 
idea of a uniform external photon field. A possible explanation is that an energetic particle 



^We find = O.OSpc, to be compared with = 1.5 x 10^"* cm (assuming A/bh = 5 x 10* M(t 



Bonnoli et al.l 120111) . Our distance val ue may challenge models in which the dissipation region distance 



is > a few pc (e.g., Ijorstad et al.ll2010l ). 
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-jQlO ^q15 ^q20 ^q25 

Frequency [Hz] 

Fig. 4. — SEDs constructed before (green open triangles), during (red filled circles), and 
after (blue open squares) the flare. The NIR, optical, UV, and X-ray data are corrected for 
Galactic extinction. The fit lines and model parameters are described in §|1] and in Table [H 
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ignition, taking place in an average- density photon region, causes the first optical fiare at 
MJD 55510. Subsequently, the blob moves away by cT5/{l + z) ^ 3.4 x 10^'' cm towards 
a region with a denser external photon field in which a doubling in the optical fiux can 
be followed by a stronger EC counterpart, as observed during the 7-ray enhanced emission 
at MJD 55517. Thereon, since the blob is moving in a region with enhanced density of 
external seed photons, the optical and 7-ray fiux variations have similar dynamic range (as 
observed at MJD 55520), until the blob leaves this denser region. Alternatively, to explain 
the whole behavior, we can invoke a modest variation of F. Subsequently, as observed in the 
post-fiare SED, the 7-ray emission decreases because of 1) the radiative cooling, and 2) the 
decrease of the external photon field due to the blob escaping the enhanced density region. 
The post-fiare SED parameters, therefore, are similar to the fiare ones, once evolution by 
radiative cooling is taken into account. During the fiare, the total power carried in the jet, 
Pjet, defined as 



Pjet = Lb + + Le + Lrad erg s" 



(3) 



(see also iGhisellini &: Celottil I2OOII ) where Lb, Lp, Lg, and Lrad are the power carried by 
the magnetic field, the cold proton 
respectively, is of the order of Pjet 



, the relativistic electrons, and the produced radiation. 



We can now discuss the absence of a harde r-when-brighter t r end in the 2-10 keV energy 
band. During the 18-months AGILE campaign. IVercellone et al.l (l2010a[) found a clear trend. 



in particular for fiuxes above (l-2)xl0~^^erg cm~^ s~^ . iDonnarumma et al.l (120101 ) show 
that a behavior similar to the N ovember 2010 one w as already present during the 2009 
December 2-3 7-ray fiare (see also iPacciani et al.ll2010l ). We can describe the harder- when- 
brighter trend in terms of a dominant contribution of the EC off the disk seed photons, 
EC(Disk), over the SSC component, probably due to an increase of the accretion rate. We 
note that an increase of and/or 7b would cause a softer- when-brighter trend, inconsistent 
with our findings. The constant X-ray photon index during the extreme 7-ray fiares in 2009 
and 2010 can be interpreted in terms of a balance of the SSC contribution with respect to 
the EC(Disk). If we assume that 7b increases significantly with respect to the 2007-2008 
ones (7b = 200-300 in 2007-2008, 7b = 700-800 in 2009-2010), we obtain both an increase 
of the EC(Disk) component (and the shift of the peak of its emission to higher frequencies), 
and a simultaneous increase of the SSC. The net result is a roughly achromatic increase of 
the X-ray emission. 
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^If considering also the contribution of relativistic protons, with (7) = 10^ , we obtain Pjet « 10''^ ergs 
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